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WPROWADZENIE 

Na wstępie rozważań należy zwrócić uwagę na wszechstronne wyko- 
rzystanie funkcji widma mocy. Analiza widmowa (spektralna) jest jedną 
z metod badania szeregów czasowych. Dostarcza ona informacji dotyczą- 
cych całkowitej wariancji procesu, omijając jednocześnie problem auto- 
korelacji. ! Szacowane widma mogą być również wykorzystane do roz- 
wiązywania niektórych zagadnień taksonomicznych. Porównanie funkcji 
spektralnej danych symulowanych i odpowiadających im danych rzeczy- 
wistych może stanowić podstawę weryfikacji modelu systemu ekono- 
micznego. Powyższe fakty zadecydowały o wyborze tematu pracy. Po- 
nieważ celem niniejszych rozważań jest zaprezentowanie wybranych 
zastosowań analizy widmowej w badaniach szeregów czasowych, artykuł 
ma charakter metodyczny. 

Część drugą i trzecią pracy poświęcono podstawowym informacjom 
odnoszącym się do funkcji widma mocy. Część czwarta zawiera zagadnie- 
nia związane z wykorzystaniem ocen funkcji widmowej do analizy sze- 
regów czasowych oraz do rozwiązywania niektórych problemów takso- 
nomicznych. Zaprezentowany w części piątej przykład empiryczny sta- 
nowi ilustrację omawianej metody analizy danych. 


1 Problem autokorelacji jest specyficzny dla szeregów czasowych i trudno go 
pominąć prowadząc badanie metodami klasycznymi. 
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FUNKCJA WIDMA MOCY 


Załóżmy, że іх, | jest stacjonarnym w szerszym sensie? procesem 
stochastycznym. Dla |Х, | istnieje możliwość skonstruowania funkcji zwa- 
nej widmową gęstością mocy. Funkcja ta przyjmuje następującą postać: 


f(a) <= (2.1) 


gdzie dF (w) reprezentuje udział czynnika о częstotliwościach z przedziału 
(о, оао) w ogólnej wariancji procesu. Wydaje się, że z praktycznego 
punktu widzenia największe znaczenie ma przedstawienie widmowej ве- 
stości mocy wynikające z twierdzenia М. Wienera—A. Chinczyna 3. Twier- 
dzenie to wyraża się zależnością: 
f(v) = > Ї ра) іе" dr = 2 IMO SIR ZĘ I ос) cos wrdr, (2.2) 
moe по по | | 


gdzie: o(r) to funkcja autokorelacji procesu (ж). Funkcję o(r) można 
przedstawić jako transformację Fouriera rzeczywistej funkcji F(w), która 
jest dystrybuantą tzw. widmowego rozkładu procesu. Jeżeli w miejsce 
funkcji autokorelacji wprowadzi się do wzoru (2.2) funkcję kowariancji, 
otrzymuje się tzw. widmo mocy procesu, które określone jest w nastę- 
pujący sposób: 
2 æ 
f* (w) = a J у(т) cos wr dr. (2.3) 
9 
Wariancja procesu wyraża się wówczas wzorem 4: 
6? = A f *(w)dw. (2.4) 


Widać więc, że widmo mocy procesu igformuje о tym, w jaki sposób 
łączna wariancja rozkłada się na drgania o poszczególnych częstotliwoś- 
ciach. W przypadku, gdy proces jest czysto losowy zarówno funkcja 
gęstości, jak i widmo mocy mają w przedziale (-П, П) wartość stałą. 


2 Mówimy, że proces stochastyczny jest stacjonarny w szerszym sensie (słabo 
stacjonarny), jeżeli jego wartość oczekiwana jest stała, a kowariancje zależą jedy- 
nie od różnicy t; — tj (nie zależą od wartości ty, tj). 

3 Рог. В. К. Oatnes, Г. Enochson: Analiza numeryczna szeregów 
czasowych, Warszawa 1972. | 

4 Wzór (14) wynika z tego, że (о) jako funkcja gęstości spełnia wa- 


runek Ў (бе) до = 1. 
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С. Granger i M. Hatanaka stwierdzają, że w ekonomicznych procesach 
stochastycznych ważniejszą rolę odgrywają wahania o dłuższych okre- 
sach, czyli widmo ma kształt szbliżony do funkcji gęstości rozkładu nor- 
malnego uciętego w wartości przeciętnej. 


ESTYMACJA FUNKCJI WIDMOWEJ 


Trwające już od końca XIX w. prace nad problemem szacowania funk- 
cji gęstości widmowej przyczyniły się do opracowania następujących me- 
tod estymacji tej funkcji 6: | 

1) „standardowej”, zwanej metodą Blackmana—Tukeya, 

2) metody transformacji Fouriera, 

3) metody filtracyjnej. 

Największe znaczenie w szacowaniu funkcji widmowej ekonomicznych 
procesów stochastycznych ma metoda standardowa. Ogólnie rzecz biorąc, 
metoda ta polega na obliczeniu ocen funkcji kowariancji, a następnie na 
wyznaczeniu jej transformaty Fouriera. N. Wiener i A. Chinczyn wyka- 
zali w swych pracach, że obliczanie widm oparte na funkcji kowariancji 
jest efektywniejsze w porównaniu z metodą bezpośredniej transformacji 
Fouriera. Zastosowana do estymacji widma mocy metoda standardowa 
wykorzystuje estymator postaci: 


1 п 


А -1 
“ер = py Слобо + 2 E Ak Ck cos ож}, (3.1) 
- k=1 

gdzie: n — liczba obserwacji w szeregu czasowym, 

II :j 

арб (j = 0,1... м). 

т 
т — punkt odcięcia funkcji kowariancji (liczba opóźnień czasowych), 
Ay — odpowiednie czynniki wagowe uśredniające widmo w paśmie czę- 


totliwości, którego środkiem jest w; (przyjmuje че hk+=0 Ша k>m). 
С. Granger i М. Hatanaka ” proponują obliczenie ocen funkcji ko- 
wariancji z wzoru: 


[5 ххк ЖЖ}. (3.2). 


$ Zob. С. W. J. Granger, М. Hatanaka: Spectral Analysis of Economic 
Time Series, Princeton 1964. 

s Dokładne. rozważania na temat tych metod zawiera praca J. Bendata, 
A. Piersola: Izmierienije i analiz słuczajnych processow, Moskwa 1972. 

7 Granger, Hatanaka: op.cit, s. 67. 
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Estymator (3.1) jest zgodnym oraz asymptotycznie nieobciążonym 
estymatorem funkcji widma mocy. % Jego wykorzystanie jest poprzedzo- 
ne doborem wag Лу. Należy zauważyć, ze ze względu na podział w pro- 
cesie estymacji podstawowego zakresu częstotliwości na (m+1) rozłą- 
cznych przedziałów, najkorzystniejsze byłoby zastosowanie wag o postaci 
filtru, którego funkcja przenoszenia jest równa: 


2m 


П 
1даој— = Sw < ој + 5— 
(3.3) 


g(w) = 
| 0 w pozostałych wypadkach. 


Brak możliwości znalezienia prostego filtru o funkcji przenoszenia 
wyrażonej wzorem (3.3) zadecydował o wykorzystaniu do estymacji wi- 
dma mocy filtrów zwanych „oknami”. 9 

Z zastosowaniem charakterystyk częstotliwościowych „okien” wiąże 
się zjawisko tzw. przecieku widma. Pojawia się w takim przypadku mo- 
żliwość przenikania części mocy związanej z częstotliwością o dużym 
znaczeniu, ale różnej од wy do oceny widma dla częstotliwości ок. W lite- 
raturze z zakresu analizy spektralnej spotyka się różne propozycje 
„okien” widmowych.10 W praktyce wykorzystuje się najczęściej nastę- 
pujące wagi: 


пк 


1) R. Blackmana—J. Tukeya Ak = 1 — 2а + 2acos (3.4) 
т 
2 [1] 11 1 ИК 
2) В. Т. „Hanning Ak = —(1 + cos —) (3.5) 
2 т 
5, ПК 
3) В. Т. „Hanning” 12 Ак = 0,54 + 0,46 cos —— (3.6) 
1 
4) Е. Parzena I NIESIE (ŻY (3.7) 


s Dowody tych własności zawarte są m.in. w pracy Т. W. Andersona: 
The Statistical Analysis of Time Series, New York 1971. 

9 Nazwa „okno widmowe” wiąże się z tym, że filtry te dają wgląd jedynie 
w pewien wycinek zbioru częstotliwości. 

10 Zagadnieniem doboru wag zajmowali się w swoich pracach m.in.: С. Jen- 
kins, E. Parzen, R. Blackman i J. Tukey. 

п Waga ta została uzyskana po podstawieniu do wzoru (3.4) а=0,25. 

12 Wagę tę otrzymuje się wstawiając do wzoru (3.4) а=0,23. 
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Należy podkreślić, że estymatory widmowe uzyskane w wyniku za- 
stosowania powyższych wag Ах charakteryzują się podobnymi własnościa- 
mi. Ciekawą systematykę charakterystyk 13 pozwalających porównywać 
efektywność wykorzystywanych okien widmowych można znaleźć w pra- 
cy A. Sokołowskiego. 14 

Szczególny wpływ na oceny widma mocy wywiera dobór punktu od- 
cięcia m (liczba opóźnień czasowych w funkcji kowariancji) oraz liczby 
obserwacji n. Dotychczasowe badania wykazują, że dobre oszacowania 
funkcji spektralnej otrzymywane są, gdy 100<n<200. Właśnie taki wnio- 
sek wysunięty został przez С. Grangera 15 i а. Jenkinsa 18. Ci sami auto- 
rzy stwierdzają, że przy wyborze liczby opóźnień czasowych powinno 
brać się pod uwagę jedną z następujących zależności: 


m <n/3, 
Т та = п/5, 
m=n/6. 


Należy również pamiętać, że w przypadku, gdy w szeregu czasowym 
występują określone wahania cykliczne liczba m powinna być wielokrot- 
nością częstotliwości odpowiadającej tym wahaniom. Wymaga się więc 
spełnienia następujących równości: 

j=(2m)/a 
gdzie: a — okres, j — liczba całkowita. 
Na zakończenie rozważań dotyczących estymacji widma mocy pro- 


cesu stochastycznego należy poruszyć problem rozkładu estymatora w 
próbie. С. Jenkins stwierdza w pracy 17, że przy założeniu normalności 


N 
rozkładu zmiennej X(t) każda składowa estymatora f* (w) ma rozkład 
w próbie określony w następujący sposób: 


ЖА 
BACH TE (3.8) 
LC) 


gdzie: s — liczba stopni swobody, х;5? — statystyka chi kwadrat o s stop- 
niach swobody. 


13 Wśród tych charakterystyk na szczególną uwagę zasługują: szerokość okna 
widmowego, wariancja estymatora widma, porównywalna liczba stopni swobody 
i porównywalna liczba niezależnych ocen. 

и А. Sokołowski: Metody badania stacjonarności jednowymiarowych ciq- 
gów losowych, Kraków 1977 (praca doktorska). 

5 Granger, Hatanaka: op.cit. 

в G. M. Jenkins, D. Watts: Spiektralnyj analiz і jego priłożenija, 
Moskwa 1972. 

9 Jenkins, Watts: op. си. 
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Założenie o normalności rozkładu procesu {Xt} powoduje, że często 
УМ 


zaleca się wykorzystywanie zmiennej log Ї* (о), która ma rozkład nor- 
malny, niezależny od częstotliwości w i słabo zależny od normalności pro- 
cesu wyjściowego. Warto zwrócić uwagę, że zależność (3.8) dotyczy jedy- 
nie estymatora widma mocy uzyskanego metodą standardową. J. Bendat 
i A. Piersol 18 wykazują bowiem, że w przypadku zastosowania np. me- 
tody transformacji Fouriera rozkład estymatora funkcji widmowej na- 
leży aproksymować rozkładem chi-kwadrat o 2 stopniach swobody. 
Oceny funkcji spektralnej można wykorzystać do obliczenia warian- 
cji empirycznej procesu Жа. W tym celu stosuje się wzór o postaci: 


1 Ko РАЧ А} 

% * (0) т-і f* (It) (3.9) 
= Еж (со) 19. 

5 [ 5 + a (wj) + 5 1. 


В. Blackman i J. Tukey !® wykazali w jednej z prac, że statystyka 
п 52/6? ma rozkład chi kwadrat o liczbie stopni swobody określonej 
w przybliżeniu jako: | 


ЖА A 
FO) m-i N rc), 
ағы ‘ a AED 2 l 
k = ——— - Xin). (3.10) 
ФР mer A (eye 
eo R [Fo + — 
2 ізі 2 


ZASTOSOWANIE OCEN WIDMA MOCY DO ANALIZY SZEREGÓW CZASOWYCH 


Zaprezentowane w pierwszych częściach pracy rozważania odnoszące 
się do podstaw analizy spektralnej oparte były na założeniu dotyczącym 
stacjonarności badanego procesu stochastycznego. Okazuje się, że szaco- 
wanie funkcji widma mocy daje dobre rezultaty zarówno w przypadku 
procesów stacjonarnych, jak i niestacjonarnych. С. Granger i Е. Parzen 
stwierdzają w swych pracach, że przy wyznaczaniu widma mocy nie jest 
konieczne badanie, czy dany szereg jest stacjonarny. Zdaniem E. Parzena, 
istnieje nawet możliwość określenia stacjonarności bądź niestacjonarności 
procesu na podstawie oszacowanego widma mocy. 

Załóżmy, że dokonano estymacji funkcji spektralnej metodą stan- 
dardową. Na podstawie wzoru (3.9) łatwo zauważyć, że analiza spektralna 


8 Bendat, Piersol: op. cit. 
18 Рог. T. H. Naylor: Modelowanie cyfrowe systemów ekonomicznych, 
Warszawa 1975, s. 355. 
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daje możliwość dekompozycji całkowitej wariancji badanego procesu na 
sumę składników odpowiadających pewnemu okresowi i przedstawionych 
za pomocą miary amplitudy. Widać więc, że strata informacji związana 
ze specyfiką szacowania widma mocy (z dużej liczby obserwacji (n) uzy- 


Ж 

skuje się jedynie (m+1) ocen estymatora f * (о), jest tylko pozorna. Po- 
jawia się bowiem możliwość określenia, które z częstotliwości wj (j= 
0, 1, ... m) w decydującym stopniu wpływają na poziom wariancji pro- 
cesu. Funkcja widmowa dostarcza również cennych informacji odnośnie 
do poszczególnych składników 2° danego szeregu czasowego. W termino- 
logii spektralnej składniki te stanowią realizacje wahań harmonicznych 
o różnych częstotliwościach. Trend jest realizacją wahań o czestotliwos- 
ciach z przedziału (0, П/2та), gdzie m jest punktem odcięcia funkcji ko- 
wariancji. Drgania harmoniczne o okresie M/k (M jest parzystą liczbą 
równych podokresów mieszczących się w okresie rocznym, a k jest nu- 
merem harmoniki (k= 1, 2, ..., M/2)) określają składnik sezonowy danego 
procesu. Składnik przypadkowy nie ma ograniczonego zakresu, a jego 
wkłady do widma mogą być rozciągnięte w całym przedziale częstotli- 
wości (0,11). Powyższe rozważania wskazują na korzyści wynikające 
z bezpośredniego zastosowania analizy widmowej do badania ekonomicz- 
nych procesów stochastycznych. Okazuje się, że oceny widma mocy mogą 
być również wykorzystywane do rozwiązywania jednego z problemów 
taksonomicznych 21, a mianowicie do grupowania obiektów jednocecho- 
wych metodą odległości widmowej. 22 


Niech $ (о) oraz ы (w) będą estymatorami funkcji widmowej uzy- 
skanymi 2 dwu różnych realizacji pewnego procesu X;. Do weryfikacji 
hipotezy: 

Ho : fx*(0)=fi* (о) 
wykorzystać można statystykę zaproponowaną przez J. Bendata i A. Pier- 
sola 23, wyrażoną następującym wzorem 24; 


3 z 2 am i (6) 
бұ S ee 521200 A (4.1) 
f у (w р 


2 Autor ma na myśli trend, wahania okresowe i składnik losowy. 

21 Taksonomia jest dziedziną zajmującą się wyodrębnianiem zbiorów jedno- 
rodnych ze względu na określone cechy. 

22 Metoda ta została zastosowana po raz pierwszy przez A. Sokołowskiego w 
pracy: Metody badania stacjonarności jednowymiarowych ciągów losowych. 

2 L. Dziembała, К. Zadora: Zastosowanie analizy widmowej do ba- 
dania wahań cyklicznych, „Przegląd Statystyczny”, 1971, nr 19. 
4 Statystyka Dy, nazwana została przez A. Sokołowskiego odległością wid- 
mową. 
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gdzie: (m +1) — liczba punktów estymacji, ny, nj — liczba stopni swobody 
odpowiednio dla 


NN 
Го) i Еко). 


Statystyka (4.1) ma rozkład chi kwadrat 25 о (m+1) stopniach swobody. 
Należy podkreślić, że jeśli obliczona wartość Пк? spełnia nierówność: 


D? < X2(m + Па, (4.2) 


(gdzie: у 2 jest wartością odczytaną z tablic rozkładu chi kwadrat dla 
(та +1) stopni swobody oraz poziomu istotności а), to przy poziomie istot- 
ności a nie ma podstaw do odrzucenia hipotezy Ho. 

Załóżmy, że wykorzystujemy do rozważań p różnych estymatorów 
funkcji widmowej. Dla każdej pary tych estymatorów wyznaczamy war- 
tości Dy? (К, 1=1, 2, ... p). Otrzymaną tą drogą macierz odległości widmo- 
wych D przekształcamy na macierz binarną P, w której O oznacza, że 
dwa estymatory są równoważne, natomiast 1, że istotnie różnią się od 
siebie. Wykorzystując w kolejnym etapie postępowania metodę elimi- 
nacji wektorów 28, dokonujemy takiego podziału analizowanego zbioru, 
aby w każdej z wyodrębnionych grup znalazły się tylko elementy parami 
do siebie podobne (w sensie kryterium odległości). Przedstawiona metoda 
może być również wykorzystana do wyodrębniania zbioru realizacji da- 
nego procesu stochastycznego. Badanie tego procesu staje się w takim 
przypadku bardziej wiarygodne (omija się weryfikację hipotezy o tzw. 
ergodyczności). 27 

Analiza widmowa zasługuje na szczególną uwagę również z innych 
względów. Stanowić ona może metodę analizy danych generowanych 
w eksperymentach symulacyjnych 28 w oparciu o modele ekonomiczne. 
Dzięki funkcji widmowej można porównywać symulowane wyniki w 
przypadku stosowania dwóch lub wielu alternatywnych polityk gospo- 
darczych. Porównanie widma danych symulowanych i odpowiadających 
im danych rzeczywistych pozwala dokonać weryfikacji modelu systemu 
ekonomicznego. 


25 Uzasadnienie faktu, ze Dy? ma rozkład chi kwadrat o (m+1) stopniach swo- 
body można znaleźć w pracy: Bendat, Piersol: op. cit. 

26 Dokładne rozważania na ten temat zawarte są w pracy Sokołowskiego: Me- 
tody badania stacjonarności jednowymiarowych ciągów losowych. 

27 Proces jest ergodyczny, jeżeli rezultaty jego badania na podstawie jednej 
realizacji można uogólnić na cały proces. 

28 Symulacja jest techniką numeryczną, służącą do opisu przy pomocy ma- 
szyny cyfrowej zachowania złożonego systemu w ciągu pewnego czasu. 


у 
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ZASTOSOWANIE ANALIZY WIDMOWEJ DO BADANIA KSZTAŁTOWANIA SIĘ CEN 
WOLNORYNKOWYCH OWSA W POLSCE W LATACH 1959—1973 


Do badania procesu kształtującego ceny wolnorynkowe owsa w Polsce 
w latach 1959—1973 wykorzystano dane pochodzące z opracowania GUS 
„Statystyka ceń”. Źródłem informacji o cenach uzyskiwanych przez rol- 
ników w transakcjach wolnorynkowych 29 jest miesięczna sprawozdaw- 
czość GUS. Przeciętne ceny dla każdego z 17 województw obliczono jako 
średnie arytmetyczne wszystkich notowań z terenu danego wojewódz- 
twa. W badaniach dysponowano 17 szeregami czasowymi, z których każ- 
dy zawierał 180 obserwacji. Ze względu na rozmiary tych- szeregów 
w opracowaniu nie przytacza się zebranych danych. 

Przed przystąpieniem do analizy kształtowania się cen wolnorynko- 
wych w Polsce w latach 1959—1973 postawiono hipotezę, że wśród 17 wo- 
jewództw występują grupy, w których dane zjawisko jest generowane 
przez ten sam proces stochastyczny. Do weryfikacji takiej hipotezy moż- 
na wykorzystać m. in. metodę odległości widmowej. W pierwszym etapie 
badania przystąpiono więc do szacowania widm mocy szeregów czaso- 
wych, reprezentujących przebieg procesu w poszczególnych wojewódz- 
twach. Estymację funkcji widmowej przeprowadzono Ша m=24, gdyż 
spodziewano się wystąpienia w analizowanych szeregach głównie wahań 
sezonowych. Do wzoru (3.1) zastosowano „оКпо" Hanninga ze względu 
na jego małą szerokość (własność ta zapewnia niską korelację między 


R (ој) a f* (w; 3) oraz małe prawdopodobieństwo przecieku). % 

Rezultaty szacowania widma mocy w poszczególnych województwach 
zamieszczono w tabeli 1. Analiza ocen funkcji spektralnej pozwala ziden- 
tyfikować cykle występujące w badanych szeregach czasowych. Prze- 
prowadzone badania nasuwają następujące wnioski: 

1) brak cyklu rocznego zauważa się jedynie w województwach w daw- 
nym obszarze: warszawskim, łódzkim, lubelskim, kieleckim i białostoc- 
kim, 

2) w przeważającej części województw występuje cykl półroczny oraz 
kwartalny, 

3) potwierdza się stwierdzenie, że w przypadku ekonomicznych sze- 
regów czasowych szczególny wkład do wariancji wnoszą wahania o dłuż- 
szych okresach. 

Należy zwrócić uwagę, że województwa, w których stwierdzono brak 
cyklu rocznego w badanych latach sąsiadowały ze sobą. | 


2 Transakcje wolnorynkowe obejmują obroty sąsiedzkie oraz sprzedaż towa- 
rów ludności miejskiej. 

о W przypadku okna „Hanninga” tzw. listki boczne nie stanowią nigdy wie- 
cej niż 2% wierzchołka głównego. 
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Kolejnym etapem badania było wyodrębnienie w zbiorze 17 woje- 
wództw grup charakteryzujących się jednorodnością (w sensie kryterium 
odległości). Wykorzystano w tym celu metodę odległości widmowej dla 
m= 24. Wyniki obliczonych tą drogą charakterystyk бұ? (К, 1=1, 2, ... 17) 
określonych wzorem (4.1) zawarto w tabeli 2. 

Ponieważ otrzymana w ten sposób macierz charakteryzuje się syme- 
trycznością, wypełniono jedynie części pod główną przekątną. Ок? ma 
rozkład chi kwadrat z (m+1) stopniami swobody. Z tabeli 31 rozkładu у? 
przy danym poziomie istotności а=0,01 32 i dla 25 stopni swobody od- 
czytano wartość krytyczną ya? wynoszącą 44,31. Następnie utworzono 
macierz binarną o 17 wierszach i 17 kolumnach (p=17), przyjmując w tej 
macierzy 0, gdy Dy? < ха?, natomiast 1 przy spełnieniu nierówności prze- 
ciwnej. W celu wyodrębnienia jednorodnych podzbiorów, w całym zbio- 
rze województw zastosowano do otrzymanej macierzy metodę eliminacji 
wektorów (metoda ta polega na kolejnej eliminacji punktów empirycz- 
nych niepodobnych do największej liczby punktów pozostałych na da- 
nym etapie rozważań). Wyniki delimitacji zawiera tabela 3. 


Tab. 3. Grupy województw wyodrębnionych metodą odległości widmowej przy 
m=24 
Groups of voivodeships selected by the spectrum distance method at m=24 


Numer grupy I II III IV У 
Митег 
województwa 1,2,5, 8, 10, 12, 17 15, 6, 7, 4 11, 16 9, 14 3, 13 


Żródło: Obliczenia własne. 


Znalezienie interpretacji dokonanego podziału wymagałoby gruntow- 
nej analizy merytorycznej przebiegu omawianego zjawiska w poszczegól- 
nych województwach. Wykracza to jednakże poza ramy niniejszych roz- 
ważań. 

Na przykładzie procesu generującego ceny wolnorynkowe owsa w Pol- 
sce w latach 1959—1973 zaprezentowano niektóre możliwości wykorzy- 
stania analizy spektralnej do badań ekonomicznych. Uzyskane rezultaty 
potwierdzają. że w przypadku gdy badacz dysponuje długimi szeregami 
czasowymi, analiza widmowa daje interesujące rezultaty. Na zakończenie 
warto podkreślić, że tematem ciekawego opracowania z omawianego po- 
wyżej zakresu może być zastosowanie funkcji widmowej do badań symu- 
lacyjnych. 


1 R. Zieliński: Tablice statystyczne, Warszawa 1972. 

33 Przyjęcie takiego poziomu istotności ogranicza możliwość znalezienia się 
w wyodrębnionej grupie elementu znacznie różniącego się od pozostałych (w sen- 
sie kryterium odległości). 
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РЕЗЮМЕ 


В статье представлены некоторые методы анализа временных рядов при 
использовании спектральной функции. Автор обращает внимание на эту соб- 
ственность спектрального анализа, которая даёт возможность определить уча- 
стие отдельных частот в дисперсии исследоваемого процесса. Много места 
посвящено также методу спектрального расстояния, который позволяет разре- 
шить одну из таксономических проблем. Представленный в работе эмпирический 
пример подтвердил пригодность упомянутых методов анализа "временных рядов. 


SUMMARY 


In the article the methods of time series analysis which are connected with 
the spectral function are presented. The author pays a special attention to this 
property of spectral analysis which lets define participation of individual fre- 
quencies in variance of the process. This article contains description of so-called 
"spectral distance method”, which is based on the comparative studies of spectra 
derived from two different realizations of the stochastic process. The methods 
described in the article were verified as exemplified by the process which generat- 
ed oat market prices in Poland in the years 1959—1973. 


